Ferromagnetism and non-local correlations in the Hubbard model 
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We study the possibility and stability of band-ferromagnetism in the single-band Hubbard model 
for the simple cubic (SC) lattice. A non-local self-energy is derived within a modified perturbation 
theory. Results for the spectral density and quasiparticle density of states are shown with special 
attention to the effects of k-dependence. The importance of non-local correlations for the fulfillment 
of the Mermin- Wagner theorem is our main result. A phase digram showing regions of ferromagnetic 
order is calculated for the three dimensional lattice. Besides, we show results for the optical con- 
ductivity and prove that already the renormalized one-loop contribution to the conductivity cancels 
the Drude peak exactly in case of a local self-energy which is not anymore true for a non-local 
self-energy. 



I. INTRODUCTION 

Band-fcrromagnctism is bound to the existence of per- 
manent magnetic moments belonging to itinerant elec- 
trons in a partially filled conduction band^. Archetyp- 
ical representatives are the classical id ferromagnets 
Fe, Co, Ni. The microscopic interpretation of band- 
ferromagnetism is one of the most fundamental and also 
most complicated many-particle problems in condensed 
matter physics. It is expected to be due to the interplay 
between ordinary, spin-independent Coulomb interaction 
(strong and strongly screened) and kinetic energy in the 
frame determined by the Pauli principle. 

A minimal model for the investigation of band- 
ferromagnetism was proposed independently by 
Hubbard^, Kanamori^ and Gutzwiller''. Despite its 
simple appearance, the Hubbard model Hamiltonian 
forms a highly non-trivial many body problem that 
cannot be treated rigorously for the general case. 

A big step forward in the understanding of correla- 
tion effects in the Hubbard model was the invention of 
the 'Dynamical Mean Field Theory' (DMFT), which be- 
comes an exact theory in the limit of infinite lattice 
dimensions^"—. The DMFT maps the lattice problem 
onto an effective single-impurity Anderson model (SIAM) 
which can be solved numerically essentially exactly by use 
of e.g. Quantum Monte Carlo methods§. 

One shortcoming of the otherwise highly successful 
DMFT is the locality (wave- vector independence) of the 
electronic self-energy strictly valid only for d — oo. So 
it may be questionable, e.g., whether such a self-energy 
is sufficient to describe angle-resolved photoemission re- 
sults. Recent efforts have therefore been focussed on re- 
gaining a certain degree of non-locality in the DMFT 
self-energy2rJ^. 

There are other approaches to the non-locality of the 
self-energy at low dimensions d = 2, 3. Coming from the 
weak coupling limit, Schweitzer and CzychoU proposed a 
method for solving the highly involved wave vector sum- 
mations that appear already in second order diagram- 
matic perturbation theoryi^. Kakehashi and Fulde used 
a projection operator method combined with the coher- 
ent potential approximation for an investigation of the 



non-local excitation spectraii^. 

Concerning ferromagnetism the few exactly known re- 
sults for the Hubbard model are of great value and can 
be used as a test frame for approximate theories. The 
Nagaoka theorem^^ states, that a saturated ferromag- 
netic order is the ground state for U = oo when one 
hole/electron is introduced into the half filled band for 
the simple cubic (SC) lattice in three dimensions (3D). 
The Mermin- Wagner theoreroii rules out ferromagnetic 
and anti-ferromagnetic order in the Hubbard model in 
dimensions d < 2 for finite temperatures p^§ - [2C| . For the 
infinite dimensional SC and FCC lattice the existence of 
ferromagnetism was proved by DMFT calculations |2l| - 

Apart from these rigorous results several works have in- 
vestigated the possibility of ferromagnetism in the Hub- 
bard model within an approximation. DMFT calcula- 
tions where done for the 3D SC and FCC lattice (2311251 
and the infiuence of next-nearest-neighbour hopping was 
investigated in [26| . Ferromagnetism in various lattices 
was investigated with a spectral density approach (SDA) 
self-energy in |27ll28| . Variational methods have been 
used in [29l - l3]| . 

A general trend can be read from these calculations. 
Two main ingredients favor ferromagnetism in the Hub- 
bard model. An asymmetric density of states (DOS) (e.g. 
the FCC DOS) and non-bipartite lattices with frustra- 
tion in the anti-ferromagnetic correlations, which can be 
generated by introducing next-nearest neighbor hopping 
t' . This shows the competitive character of ferro- and 
anti-ferromagnetic correlations in the Hubbard model. 

In this paper we investigate the infiuence of non-local 
correlations on ferromagnetic order in the SC lattice. To 
this end we shall apply the Modified Perturbation Theory 
(MPT), which was originally used only for solving the 
SIAM within the DMFT procedure, directly to the full 
Hubbard problem. The MPT leads to an explicitly wave- 
vector dependent self-energy which decisively determines 
the single-electron spectral density. It is well known that 
the latter provides the bare line shape of a spin and angle- 
resolved (direct or inverse) photoemission experiment. 

The paper is organized as follows. In section [H] we 
introduce the Hubbard model Hamiltonian, derive the 



MPT self-energy and discuss its properties. Then ther- 
modynamic quantities as the paramagnetic static suscep- 
tibihty and the optical conductivity are derived. 
In section mil the numerical methods for dealing with the 
complicated momentum summations are presented. 
Section IIVI contains the results and interpretation of our 
numerical calculations. Finally we give a summary and 
conclusion in section IVl 



II. THEORY 



mogeneous external magnetic field {z^i = ±1). We have 
chosen the band center of gravity to = and the hopping 
strength t such that the free electronic bandwidth W is 
equal to one throughout the paper. An (approximate) 
solution of the Hubbard model is found, if we are able to 
calculate the electronic Green's function (GF): 



Gka(S) - (S + M - e(k) - I]k.(£;)) 
or more precisely the electronic self-energy 'S-^criE). 



(2) 



Hubbard model 



B. self-energy 



The Hamiltonian of the Hubbard model is given by: 
H = t^ cl^Cja+^{zaB+to)nicr + U^ni^fni^l. (1) 

{Lj)a ia i 

Here t denotes the nearest neighbor hopping strength 
(the sum over Rj extends only over the nearest neighbors 
of Ri), U is the local coulomb repulsion and B an ho- 



lt is now a well-known fact, that the self-energy of the 
Hubbard model becomes a purely local quantity in the 
limit of infinite dimensions (d — ?> 00)^1^. However the k- 
dependence will certainly play a crucial role in the more 
realistic case of d = 2,3. In the weak coupling limit 
{U <^ W) the second order perturbation theory (SOFT) 
is a good starting point for the investigation of non local 
correlation effects. The SOFT self-energy is given by^-: 



J 



Sk,(i?) = 4^^) + s. 



{SOC) 



(E) 



(3) 



U{n^,) + U^Y1 ^*'' / '^^SnAx) J dySn-Ay) J dz5_R_,(z) 

Ft, 



f-{x)f-{y)f-i-^) + f-{-x)f-{-y)f-i^) 

E + iO+-x-y + z 
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The sum extends over all lattice sites R. Schweitzer and 
CzychollJA gave a method for the calculation of ([3]) by 
collecting all symmetry equivalent points in shells and 
recast the sum over lattice sites into a sum over shells. 
They showed, that this sum can be truncated after a 
finite number of shells. However, their method of calcu- 
lating the real and imaginary part of ([3]) is still numer- 
ically very demanding. We will show in section IIIII how 
to speed up the computation to allow fully self consistent 
calculations for arbitrary band fillings n. There is a cer- 
tain arbitrariness in ^ concerning the spectral densities 
(SD) appearing in the formula. In an expansion strictly 
to order U^ the free SD has to be chosen^. But one 
could also renormalize the theory by using the full SD 
in a self-consistent manner. It turns out, that only the 
first choice will reproduce certain exact results^^. To be 
specific we give the form of the SD used: 
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(4) 



Ha is fixed by the condition, that the free occupation 
number is equal to the full occupation: (n^)'"' = (na)- 

Notice, that this choice of /ijj is equivalent to a SOFT 
"around Hartree-Fock (HF)" at half fiUing where the full 
/i from © is taken in the HF-SD. Only this choice will 



result in the now widely accepted three peak structure 
of the density of states (DOS) and will give a "smooth" 
change of the DOS away from half filling. For more qual- 
itative discussions we refer the reader to the results sec- 
tion. 

To extend the validity of the self-energy to larger values 
of U we use the following Ansatz for a modified pertur- 
bation theory (MFT): 



Sk.(£^) = U{ny + 



«k.4f^^(i?) 



1 - 6kcrS 



(SOC) 



(E) 



(5) 



This form for the self-energy was proposed by Kajueter 
and Kotliar'^^ for the Anderson impurity model (SIAM). 
They used the first two spectral moments and an ad- 
ditional condition for the chemical potential to fix the 
parameters Oko- and foko-- This method of fixing the pa- 
rameters was afterwards modified by Fotthoff, Wegner 
and Nolting-- for the same model in order to reproduce 
the first four moments of the SD correctly. We will fol- 
low this latter approach but now applied to the full lattice 
Hamiltonian (jlj. 
To fix the appearing constants in ([SJ we use the high en- 



ergy expansion of the self-energy for the Hubbard modeh 



E" 



The first three coefficients can be obtained from the first 
four moments M^' of the SD^" via the high energy ex- 
pansion of the electronic GF ([2]) : 



GUE) = ^ E 



1 ^ M^r^ 



T71 — 



(7) 



and are given in the appendix (JA1|) . By expanding also 
the rhs of ^ we can determine the coefficients akcr,&ko- 
and get finally the MPT self-energy: 



Sk.(i?) = c/(n-.)+ (4",'^^^(£;))-i 



D 
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a 
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kCT 



(C, 



kCT ' 



(8) 

where D^^ denotes the third moment of ([3]) as given in 
the appendix (jB2|) . 

The MPT self-energy can be proved to be exact in a 
variety of limiting cases. It trivially fulfills the limits of 
[7 = and n = 0,n — 2. More interesting is case of zero 
bandwidth limit t — s- 0. A straightforward calculation 
yields: 






(E) = U{n^^) + 



U^n^,)il-{n.„)) 



E 



rB 



M 



t/(l-(n_.))' 
(9) 
which is indeed the correct form of the "atomic limit" ■^°. 
Expanding the self-energy for small U reproduces the re- 
sult of perturbation theory ([3]) with corrections only of 
order U^ . Therefore the MPT should be correct for small 
U and show e.g. Fermi liquid behaviorJ^. 
Since the first four spectral moments are reproduced cor- 
rectly by construction, they should be correct at large en- 
ergies \E\ ^ 1, too. In particular the position and shape 
of the upper (lower) Hubbard band for n < 1 [n > 1) 
will become exact in the strong coupling limit (C/ 3> 1) 
which is known to be a weak point of SOPT alone. In 
this respect the MPT self-energy is in accordance to the 
t/U strong coupling expansion of Harris and Lang o^^i^^ . 
To conclude this discussion we summarize our findings. 
We have proposed a fully k-dependent MPT self-energy 
which fulfills the atomic limit and shows reasonable be- 
havior in the weak and strong coupling region. Therefore 
there is well-founded hope that our theory will give rea- 
sonable results also in intermediate coupling region. 



C. thermodynamics and transport 

1. paramagnetic static susceptibility 

The developed theory allows for a self-consistent cal- 
culation of the magnetization in a possible appearing fer- 



romagnetic region. To test the system regarding a ferro- 
magnetic phase transition, we will calculate the param- 
agnetic static susceptibility, which is defined as follows: 



X 



'M 



(T) ^J2^Biz^{na))\T,B=0,{n^) = {n^)- 



(10) 



The zero crossings of the inverse of pHt indicate the 
points where the paramagnetic phase become suscepti- 
ble to a ferromagnetic phase transition. 

For an evaluation of (jlOp one has to perform the deriva- 
tive analytically and get after a lengthy calculation an 
explicit form for the susceptibility as a functional of the 
(self-consistently determined) paramagnetic self-energy. 
Since the expressions are rather long, we do not give them 
here. 



2. optical conductivity 

The optical conductivity in linear response is given by 
the retarded current-polarization GF— : 



Ba 



{E) = -{{f-P'^ 



(11) 



where a, j3 denote the Cartesian coordinates of the opera- 
tors. By writing down the EQM of this GF and exploiting 
the connection j 



i-^[P, if]_ this can be rewritten as: 



^/3a(£;) = _iL?!ii!!Ll + ,vvi^^i^. (12) 



E 



E 



For a tight binding (nearest neighbor hopping) model 
the operators are given as P = g ^^ ^ RifT-io- and j = 
~lv^I^(im),CT(^«^^™)Ci!TCmCT- With these operators the 



first term of (|12l) can be calculated and in case of a simple 
cubic lattice simplified to give the zero frequency Drude 
weight of conductivity: 

Re(a^"(i? + zO+)) = -7r^a^<5(i?)^ ^cos (fc„)(nk.). 

kcr 

(13) 
The second term in p^ is the current-current GF. It 
represents the influence of electronic correlations. We 
approximate this GF on the 'one loop' level in a dia- 
grammatic expansion. The explicit calculation is given 
in appendix ([Cj and it is shown that the real part con- 
sists of two parts. One is proportional to 5{E) and can- 
cels the Drude peak in case of a local self-energy exactly. 
The second term yields: 



Re(cr^"(£; + iO+)) = ^ap-rj- E WkaWk/3 

kCT 



/ 



dxSUE + x)SUx)l^^^^^J^-^^±^, (14) 



with Uka = 9fe„ek- 

Note, that this approximation becomes a rigorous re- 
sult in infinite dimensions because the self-energy is a 



local quantity and all higher vertex corrections vanish in 
this case^. 

We will only show results for the contribution from 
(fT4|) hoping that the exact cancellation is retained in the 
k-dependent case approximately at least. The neglect 
of vertex corrections can not be justified rigorously in 
case of a k-dependent self-energy, because Ward iden- 
tities may be violated. We show therefore only results 
for the 3D SC lattice, where the k-dependence of the 
self-energy is not so pronounced (especially at the Fermi 
level) and refer the reader to the more specialized litera- 
ture for a thorough discussion of this point^. 



where F{x,y, z) denotes the product of Fermi functions 
in ([3]). This twofold convolution can be solved very ef- 
ficiently by a fast Fourier transform. From this Y.^^ (E) 
can be obtained in the complex plane via the spectral rep- 
resentation of Greens functions^. The free shell density 
of states needed for the calculation we have computed 
and stored beforehand to very high precision. The num- 
ber of shells necessary to get converged results depends 
strongly on the coordination number of the underlying 
lattice. After some tests we have used 201 shells for the 
2D and 61 shells for the 3D lattice for all calculations in 
this work, which was sufficient to get mostly converged 
self-energies. 



III. COMPUTATIONAL METHODS 

The self-energy (j8]) is a functional of the chemical po- 
tential and different correlation functions. It has to be 
calculated self-consistently. To this end we need a fast 
way of calculating integrals of the form: 



(...) 



-E 



d£;/_(^)Fk(^)5k(^), (15) 



where f-{E) is the Fermi function, F\^{E) a polynomial 
of low order in E and Su{E) the full SD. It is hopeless 
to perform the four dimensional integral directly because 
the SD is a strongly peaked function. However one can 
replace the energy integration by a sum over the poles of 
/_ {E) in the upper complex plane. The usual Matsubara 
form of the Laurent expansion of /_ (E) is not suitable 
here, because it converges very slowly with an increasing 
number of poles. Recently Ozaki 40] proposed a different 
pole expansion for the Fermi function which gives a good 
approximation for a large energy domain down to very 
low temperatures with only a few hundred poles. We use 
this expansion for a numerical very accurate determina- 
tion of the energy integral. The remaining k integration 
over the irreducible wedge of the Brillouin zone can then 
be performed directly. 

For the exposed procedure we need the values of the 
SOFT self-energy at the Ozaki poles in the complex 
plane. For their determination we first rewrite (|31) as 
a sum over shells of symmetry equivalent points: 



^SOC) 



{E)^U 



N 

1 Y^ 






GC>yS^\e), 



(16) 



^{o) 



where G\^ denotes the k dependent shellfactor given by 
the sum over the exponentials in (|3]) within one shell 
and Y^a (E) the remaining energy dependent part. The 
imaginary part of the latter is given by: 

Iml](,")(£;) = -TT I dx I dy (17) 

Sil\x)si%{y)si%{x + y - E)F{x, y,x + y~E), 



IV. RESULTS AND DISCUSSION 

Our k-dependent self-energy allows the description of 
homogeneous phases (paramagnetic/ferromagnetic). We 
will discuss the electronic properties in the paramag- 
netic state of the three dimensional Hubbard model first, 
knowing well that in certain parameter regimes (e.g. near 
half filling at low temperatures) anti-ferromagnetism is 
expected in principle. This restriction is shared with 
other self-energy approaches (e.g. DMFT calculations) 
and a comparison to these should render our results use- 
ful. 



A. QDOS and spectral density 

In Fig. [1] we show the quasiparticle density of states 
(QDOS) at low temperature (T = lOK) for two dif- 
ferent band fillings and various interaction parameters 
U. The three peak structure of the QDOS is clearly 
visible. Upper and lower Hubbard band are roughly 
separated by the coulomb interaction strength U and 
there appears a Kondo resonance near the Fermi level 
{E = eV in figures). By increasing U, this resonance 
decreases but stays finite also for large U. As a conse- 
quence there is no clear metal-insulator transition (MIT) 
for n = 1 as is found in DMFT calculations. To il- 
lustrate this further we show the inverse effective mass: 
^ = iw Ek(l -^Sk(0))~' in Fig.H Although there is 
no clear transition point the system should be insulating 
above U/W w 10 due to the large effective mass of the 
quasiparticles. This finding is in qualitative agreement 
with the non-local theory of Kakehashi and Fuldei^. 
The comparison of the local approximation (only zeroth 
shell of the SOC is taken into account) and the full k- 
dependent calculation shows decisive effects of the latter. 
Whereas the local theory fulfills the Luttinger theorem 
(the QDOS at the Fermi level is equal to the free DOS), 
the k-dependence leads to a reduction of states at this 
point. This is understandable because the Luttinger- 
Ward argument only holds for local self-energies^i^. 
Another effect of the fc-dependence are the peaks in the 
upper (n = 0.75, 1.0) and lower (n = 1.0) Hubbard band 
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FIG. 1: (color online) QDOS for various interaction parame- 
ters U/W; full lines: k-dependent self-energy results; broken 
lines: local self-energy and interaction free result; Parameters: 
T = lOK, upper panekn — 0.75, lower panehn — 1.0 




FIG. 2: Inverse effective mass as function of interaction pa- 
rameter U/W; Parameters: T = lOK, n = 1.0 



for intermediate coupling {U/W « 1). For interpreta- 
tion of these peaks we have plotted the spectral den- 
sity together with the imaginary part of the self-energy 
along special directions within the first Brillouin zone 
for n = 1 in Fig. |31 The self-energy shows typical 
Fermi liquid behavior. The imaginary part is near zero 




FIG. 3: Spectral density (upper figure) and imaginary part 
of MPT self-energy (lower figure) along special directions 
in the first Brillouin zone of the simple cubic 3D lattice: 
M{tt, TV, 0) -^ X{tv, 0, 0) -^ r(0, 0, 0) -^ R{-K, TT, tt) -s- M; Pa- 
rameters: n = 1, T = lOK, U/W = 1 



(zero only at T = 0) at the Fermi level and decreases 
quadratically with increasing energy. This leads to in- 
creasing damping effects in the dispersion of the Kondo 
resonance particularly strong at the F and R points. 
At these points the self-energy shows a strong enhance- 
ment around E = — / -|- 0.5 respectively and increases 
then abruptly to zero for lower/higher energies. There- 
fore we find no damping effects in this energy region 
and quasiparticle states with energies lower/higher than 
the threshold energy will have infinite lifetimes. This 
is clearly visible in the spectral density where in the 
lower/upper Hubbard band "bridges" of sharply peaked 
states appear at both k points. Similar effects where 
found in SOFT calculations by Schweitzer and CzychoU— 
and in the projection operator method of Kakehashi and 
Fulde^^ but here they are much more pronounced. The 
reason is most likely the increased numerical accuracy of 
our calculation (see section llll|) . 

Comparing our results to a recent dynamical cluster ap- 
proximation^^ (a cluster extension of the DMFT in or- 
der to retain a k-dependence) there is good agreement 
below the MIT {U/W = 0.67 in the mentioned work). 
The dispersion of the Kondo resonance and the maxima 
of the lower/upper Hubbard bands show essentially the 
same behavior as in our calculation. Above the MIT 
{U/W = 1) there are of course discrepancies because of 
the missing MIT in the MPT. 

With increasing temperature damping effects will become 
more important. This is shown in Fig.|4l The Kondo res- 
onance peak is diminished with raising temperature and 
tend to vanish completely at higher temperatures. The 
reason for this can be found in the inset of Fig. 01 With 
increasing temperature the averaged imaginary part of 
the self-energy at the Fermi level decreases starting from 
zero at T = QK . For low temperatures a typical Fermi 
liquid behavior is obtained (^ T^). 
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FIG. 4: (color online) QDOS (Kondo resonance) at various 
temperatures T\ Inset: averaged imaginary part of MPT self- 
energy at Fermi level E = Q\ Parameters: n= \, U/W — 2. 



FIG. 5: Spectral density (upper figure) and imaginary part of 
MPT self-energy (lower figure) along special directions in the 
first Brillouin zone of the simple cubic 2D lattice: M{t:, n) — >■ 
r(0,0) -^ X{tt,0) -s- M; Parameters: n = 1, T = lOK, 
U/W = l 



For the two dimensional SC lattice the effects of a non- 
local self-energy should be more drastic then in 3D as a 
direct consequence of the reduced coordination number. 
Fig. [5] shows the spectral density and imaginary part of 
the self-energy of the SC 2D lattice at half filling. We find 
again states with infinite lifetime at F and Af point due 
to the vanishing imaginary part of the self-energy. The 
self-energy shows Fermi liquid behavior ('^ £^^) at large 
portions of the Brillouin zone. However, this behavior 
is changed near the Fermi surface. At the X point and 
the midpoint between F and M (f , f ) the self-energy 
shows a linear energy dependence. This "marginal Fermi 
liquid" behavior is a direct consequence of the perfect 
nesting properties (e(k) w e(k + Q), where Q is a re- 
ciprocal lattice vector) of the 2D Fermi surface^"^ at 
half filling. Electrons can be scattered efficiently due to 
the large phase space. This leads to a second effect in 
the spectral density - the formation of "shadow" bands. 
They are visible as dark lines running from the special 
points X, (fjf) to F, M with a slope determined by 
the condition E « e(k-f Q). These shadow features, 
which have already been described by Vilk^, are no real 
quasi particle band but merely thermal excitations cor- 
responding to a local minimum of the imaginary part of 
the self-energy. To illustrate this point we have plotted 
the spectral density (upper panel) and self-energy (lower 
panel) at k = (||, |^) in Fig. El The green dash-dotted 
line obeys E + fi — e(k) and its crossing points with the 
real part of the SE (black line) define the quasiparticle 
excitations of the system. In the SD only one of the 
three excitations forms a peak (near E = 0), the others 
are strongly damped by a large imaginary part of the SE 
(red line). The black dashed line marks the position of 
the shadow band. There is no real excitation energy but 
we find a local minimum of the imaginary part of the SE 
which leads to observed shadow feature. 
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FIG. 6: (color online) Spectral density (upper figure) and real 
part (lower figure: black line) and imaginary part (red line) of 
the self-energy for the simple cubic 2D lattice at k = (ff) |f )• 
Vertical broken line: position of the shadow band. The cross- 
ing point the green dash-dotted line with the real part of 
the self-energy marks the position of quasiparticle excitations. 
Parameters: n = 1, U jW = 1, T = lOK. 



B. conductivity 

From P^ it becomes clear, that the optical conduc- 
tivity is mainly determined by the number of available 
quasiparticle states. In Fig. [7] the inverse static conduc- 
tivity (resistivity) at half filling is shown as a function 
of temperature for U = 2. At low temperatures the re- 
sistivity increases quadratically. This results from the 
reduction of the QDOS at Fermi level as shown in Fig. S) 
The resistivity rises until the thermal energy is sufficient 
to excite electrons from the lower to the upper Hubbard 
band (fcsT « U). Then it will decrease going through 
a minimum and rise again. The inset shows the opti- 
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FIG. 7: Correlation part of the resistivity (inverse of the static 
conductivity) as a function of temperature T for the 3D simple 
cubic lattice. Inset: optical conductivity as function of E. 
Parameters: n = 1, [/ = 2; Inset: T = lOOK. 



cal conductivity at fixed temperature T — lOOK. The 
conductivity decreases with increasing energy due to the 
lack of states between the Kondo resonance and the up- 
per Hubbard band. As soon as the energy is sufficient to 
excite electrons from the Fermi level to the upper band, 
the conductivity will rise strongly reaching a maximum 
soon and decrease again. These findings for the conduc- 
tivity are in well agreement with DMFT results for U 
below the metal-insulator transition^ 



C. inverse paramagnetic susceptibility and 
ferromagnetic phase transition 

The inverse paramagnetic static susceptibility (IPS) 
(TTOt can be used as a tool for finding borders of a fer- 
romagnetic phase transition in the n-U diagram. Its zero 
crossings will mark the critical points. In Fig. [8] we show 
the IPS for the 2D and 3D SC lattice at low tempera- 
ture (T = lOK) for various U. Whereas we find zero 
crossings above a critical U/W « 15 in the 3D case, 
there is no point of phase transition in the 2D case when 
we use the full non local MPT self-energy. Only when 
using the truncated local version of the MPT we get a 
phase transition in 2D also. This shows first of all, that 
the predictions of the non local theory are in accordance 
with the Mermin Wagner theorem and secondly, that the 
non-locality of the self-energy is crucial in oder to get the 
result. We would like to mention that the IPS curve for 
U/W = 500 is saturated in the sense, that increasing 
U/W further will not change its shape drastically. 
We come now to the discussion of the magnetic properties 
of the 3D system. The rather high critical U rcficcts the 
fact, that the SC lattice is not particularly susceptible 
to ferromagnetism (the competing anti-ferromagnetism 
is not suppressed by frustration, like in the FCC lattice). 
The IPS crosses the zero axis at two points. These mark 
the lower and upper bound of the ferromagnetic region. 




FIG. 8: Inverse paramagnetic susceptibility as a function of 
band filling n for the 2D (upper figure) and 3D (lower figure) 
system calculated with the full k-dependent self-energy. For 
the 2D system the result of local self-energy is shown also 
(broken lines). Parameters: T — lOK. 




FIG. 9: Curie temperature Tc as function of band filling n 
for various U/W . Inset: relative magnetization (n^ — ^i)/^ 
as function of interaction strength U/W . Parameters inset: 
n = 0.9, T = lOK. 



The Curie temperatures for different U/W are plotted 
in Fig. |9l The maximal Tc is reached for a band filling 
of 71 « 0.91. Starting from Ucrit/W w 15, Tc increases 
quickly with increasing U/W and running into satura- 
tion for larger values of U/W. The same is true for the 
magnetization m — n-f — n^ and the phase border. The 
magnetization at n = 0.9 and T — lOK is shown in the 
inset of Fig. [9l The electron system is far from satura- 
tion, the polarization reaches ~ 16% for U/W = 70 and 
increases slowly for stronger interaction parameters. 

The full phase diagram is shown in Fig. 1101 The lower 
phase boundary is decreasing with increasing U/W to 
lower n up to U/W ~ 350 where it takes the value 
Ticrit ~ 0.845. Increasing U/W further will not increase 
the ferromagnetic region but we observe a slight shifting 
to higher n again. The upper phase boundary growths 
monotonically with increasing U/W reaching a value of 
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FIG. 10: Ferromagnetic phase diagram of the 3D simple cubic 
lattice as function of U/W and band filling n. Parameters: 
T = lOK. 



n sa 0.969 for U/W — 500. This result is in accordance 
with Nagaokas theorem. 

As expected, in the local approximation the critical U: 
Ucrit/W « 5 at T = 10 K for the 3D SC lattice is lower 
than in the full k-dependent case and agrees well with 
findings from DMFT calculations {U/W > 3)^. 



V. SUMMARY AND CONCLUSIONS 

In this work we have derived a non-local self-energy 
for the Hubbard model within a modified perturbation 
theory approach.lt was shown, that this self-energy 
fulfills a variety of limiting cases (e.g. weak coupling and 
atomic limit) and shows the correct high energy behavior 
by construction. Numerical tools for the evaluation of 
this self-energy were introduced, in particular to solve 
the complicated momentum integrations. 
We show results for the two- and three-dimensional 
simple cubic lattice, discussing in detail the influence 
of non-locality, temperature and interaction strength 
on the self-energy, spectral density and quasi-particle 
density of states. These results are then used for the 
interpretation of the calculated optical conductivity and 
resistivity curves. 

The inverse paramagnetic susceptibility was calculated, 
showing that there is no ferromagnetic phase transition 



in the two-dimensional but for the three-dimensional 
lattice. The ferromagnetic/paramagnetic phase diagram 
for the three-dimensional lattice is then constructed. 
Our findings emphasize the importance of non-local 
correlations in the Hubbard model in low dimensions, in 
particular for the fulfillment of Mermin- Wagner theorem. 

The strength of the present approach is that it allows 
self-consistent calculations at arbitrary band-fillings 
with a fully k-dependent self-energy. This comes at the 
expense of neglecting certain correlation effects. Most 
severe in this respect is the missing metal-insulator 
transition (MIT) at half-filling which should also occur 
in finite dimensions at roughly U/W « 1 as indicated 
by DMFT calculations^"^. It is interesting to note, that 
the (local) MPT, when used as an impurity solver in a 
DMFT calculation for the infinite dimensional Hubbard 
model does show a MIT in the correct U/W region, 
whereas it does not when used as a self-energy for the 
lattice Hamiltoniani^. 

Another shortcoming of the current state of the the- 
ory is the inability to check for antiferromagnetic phases, 
which ultimately should appear near half-filling. The rea- 
son for this is in the Ansatz of the self-energy ([8]), which 
only allows the calculation of homogeneous phases. This 
does not mean, that there are no possible antiferromag- 
netic solutions within the MPT approach. It would be 
an interesting task for a forthcoming work, to extend the 
MPT in this direction. An indirect hint for the presence 
of antiferromagnetic correlations within the MPT can be 
derived from the breakdown of ferromagnetic order near 
half filling, which could have its origin in competing an- 
tiferromagnetic correlations. 



Appendix A: coefficients of self-energy expansion 
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Appendix C: Derivation of the optical conductivity 

The density-density GF in (renomaUzed: free propagators are replaced by full ones) diagrammatic one-loop expan- 
sion is given by: 



^("-kcr; "-k'<T))_E„ ~ '^kk'^ ^ Gka{iErn)Gl^a{iEn + iE„i) 



P 



<^kk' / / dxdySkaix)S]^cr{y)- 
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lEn +x -y 



(5kk' / dxf.{x)Skaix) (Gk.(a; + E + iO+) + G^aix -E- iO+)) . 
From this we get for the correlation part of conductivity ([T2|) : 

Reaf;(£; + ^0+) = Re i —^^-—1 Y.^du,e,,){dk^t^) j dxf-{x)S^^{x) (Gk.(x + E + zO+) + G^^{x -E- iO+)) 

= -2qH{E)^ Y.^dk,e^){du^e^) j dxf_{x) (ImGk.(x)ReGk.(x)) 

kcr 

5^(9fe,ek)(9fe„ek) ( dxf-{x)S],.{x) (Si^^x + E) ~ 5k.(a; - E)) . 
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The first term of tliis result can be recast into tlie form of the Drude contribution (jl3|) but with opposite sign: 
-2qH{E)^Y.('^k,eu){dk^e^) f dxf.{x) (IniGk,(x)ReGk.(x)) 

kcr 

kcr " 

= ^TTqH{E) 1 ^(afc„ dk, eic) y da:/_ (a;)5k. (a;) 

kcr 

kCT 

where the second and third step is only allowed when the self-energy does not depend on k and several steps require 
a diagonal mass tensor (as it is for the SC lattice). 
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